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Abstract 

Quantum algorithms for quantum chemistry map the dynamics of electrons in a 
molecule to the dynamics of a coupled spin system. In order to reach chemical accu¬ 
racy for interesting molecules, a large number of quantum gates must be applied which 
implies the need for quantum error correction and fault-tolerant quantum computa¬ 
tion. Arbitrary fault-tolerant operations can be constructed from a small, universal 
set of fault-tolerant operations by gate compilation. Quantum chemistry algorithms 
are compiled by decomposing the dynamics of the coupled spin-system using a Trotter 
formula, synthesizing the decomposed dynamics using Clifford operations and single¬ 
qubit rotations, and finally approximating the single-qubit rotations by a sequence of 
fault-tolerant single-qubit gates. Certain fault-tolerant gates rely on the preparation 
of specific single-qubit states referred to as magic states. As a result, gate compilation 
and magic state distillation are critical for solving quantum chemistry problems on a 
quantum computer. We review recent progress that has improved the efficiency of gate 
compilation and magic state distillation by orders of magnitude. 
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INTRODUCTION 


Quantum simulation is the mapping of one quantum system to a second, well-controlled 
quantum system in order to both understand how the system changes as parameters are 
varied and to allow measurements that may be hard to access in the original system m- 
Quantum simulation can offer an advantage over classical computer simulation, since there 
is no compact way to represent arbitrary quantum states with classical bits. The application 
of quantum simulation to quantum chemistry is seen as a promising application for a future 
quantum computer. 

Calculations have shown that a full-configuration interaction (FCI) calculation on a quan¬ 
tum computer can be solved with a number of operations that is polynomial in the number 
of spin-orbitals assuming that an approximation to the quantum state of interest can be 
prepared [3j. The bulk of the algorithm consists of simulating the dynamics of the second- 
quantized molecular Hamiltonian with N spin-orbitals: 

Ff ^ ^ hjkCjCk T ^ ' hj kirn, QQQCm (1) 

j,k j,k,l,m 

where hjk are the single-electron energies corresponding to the kinetic energy and the attrac¬ 
tion to the nuclei, hjkim are the two-electron energies due to electron-electron interaction, 
and c'j(cj) is an electron creation (annihilation) operator for the spin-orbital j. The dynamics 
are simulated for a set of times t, 2 1, 4 1,... 2 k t = T where t sets the maximum energy that 
can be observed and T sets the precision due to the standard limits of the Fourier transform 

EE- 

The quantum simulation is performed by mapping the dynamics of the molecular Hamil¬ 
tonian to a series of quantum gates (Figure [l]). The first step is to map fermionic operators 
to spin operators [6]. There are two generic methods of performing this mapping: Jordan- 
Wigner, which encodes orbital occupation in individual spins at the cost of iV-body spin 
operators to account for the exchange symmetry EEli, and Bravyi-Kitaev, which uses 
log(iV)-body spin operators at the cost of losing the direct mapping from orbitals to qubits 
0 ED]. For restricted Hamiltonians, e.g., models where only spatially local interactions 
between electrons are allowed, simpler mappings are possible |6]. 

A universal quantum computer can generate an arbitrarily accurate approximation to any 
unitary evolution on 2 n two-level systems or qubits. A key result of quantum computation is 
that this unitary operation can be broken into arbitrary single-qubit gates and any entangling 
two-qubit gate mm- In this framework, the compilation of various unitary operations can 
remain challenging. 

For quantum chemistry, the simulation of each M -body spin operator is broken down 
into operations on one and two qubits, which can be combined using the Trotter formula, 
exp j- ~ (n.exp jr'Hk/m^ . Based on estimates of the required Trotter 

time-steps and the choice of spin-transformation, the total number of operations scales as 
N 8 |E3], for arbitrary h values, to N 6 using Jordan-Wigner jb|j or N 5 log(AT) using Bravyi- 
Kitaev [15] for values of h commonly seen in molecules. This is quite promising for FCI since 
the scaling within the regime of worst-case coupled-cluster singles and doubles (CCSD), IV 6 , 
[16] and the number of required quantum bits scales only as N. We also note that there is a 
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Figure 1: Presented are 
the steps to simulate the 
dynamics of a Hamilto¬ 
nian on a quantum com¬ 
puter. This review fo¬ 
cuses on recent progress 
in the final step consist¬ 
ing of gate compilation 
and magic-state distilla¬ 
tion. 


C 


Dynamics of the Second-Quantized 
Molecular Hamiltonian 


Mapping of Spin-Orbitals to Qubits 




c 

c 


Exponentiated Weighted Sum 
of Pauli Operators 


Trotterization 




Product of Exponentiated 
Pauli Operators 


Unitary Operator Decomposition 




C 


Clifford Gates and Arbitrary 
Single Qubit Rotations 


Fault-Tolerant Procedures 




C 


Gate Compilation and 
Magic-State Distillation 


A 


connection between electron correlation and the required number of Trotter-steps resulting 
in quantum algorithms with only N 4 total operations for molecules well-described by the 
Hartree-Fock approximation [ IT] . 

These estimates assume perfect qubits and operations. Actual qubits have errors requir¬ 
ing fault-tolerant methods for robust quantum computation. Fault-tolerant quantum error 
correction encodes a single logical qubit into multiple physical qubits. This allows us to 
preserve information in the logical level by detecting and correcting errors at the physical 
level H7j. Operations on logical qubits must be fault-tolerantly performed, which requires 
the careful design of logical gates from physical gates. Fortunately, only a restricted set of 
gates is required to achieve universal quantum computation [18]. For many error-correcting 
codes, the fault-tolerant gates fall into two classes: gates that can be directly applied to the 
logical qubit and gates that require the consumption of a specially prepared ancillary state, 
known as a magic state [IPj. 

The operators that were decomposed assuming no quantum error correction now need to 
be decomposed into this restricted gate set. The effect of this decomposition on the total 
number of qubits and gates has been considered previously for simulations of spin and molec¬ 
ular Hamiltonians [20H22] . but in the last few years there has been a remarkable reduction 
in the required resources. In this review, we will examine advances in the compilation of 
arbitrary single-qubit unitary gates from restricted gate sets and improvements in the gen¬ 
eration of magic states by a process called distillation, where many low-quality magic states 
are converted into a few high-quality magic states in a fault-tolerant manner. 


Quantum Gates and Teleportation 

We start with some useful notation and definitions. For further details, we recommend 
the textbook Ref. [T8]. The Pauli matrices on a single-qubit j are represented by Xj = cfj, 
Yj = dj, and Zj = crj. The two common bases are the computational or z basis corresponding 
to eigenstates of Z, Z |0) = |0) and Z 11) = — 11), and the x basis corresponding to eigenstates 
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of A", X |+) = |+) and X |—) = — | —}. The Hadamard matrix is a self-adjoint matrix that 
transforms between these bases, HXH = Z. The standard two-qubit gate is the CNOT (i,j) 
which adds the value of qubit i in the computational basis to qubit j modulo 2 and can be 
written as 

CNOT(i,i) = |0} (0|| Ij + |1) <1| ( X,- (2) 

where Ij is the identity on the jth qubit. 

Two important groups of unitary operations are the Pauli group and the Clifford group. 
The Pauli group, V, is generated by multiplication of individual Pauli matrices on each qubit. 
The Clifford group, C , is the group of unitary transformations that maps the Pauli group 
to itself; for g e C, gVg ' = V. On a single-qubit, the Clifford transformations correspond 
to the symmetry elements of the chiral octahedron (the point group O). The Clifford group 
can be generated by the single-qubit operators H, S = |0) (0| + i |1) (11, and the two-qubit 
operator CNOT. 

Each element of the second-quantized molecular Hamiltonian can be mapped to a sum of 
elements in V [23]. The dynamics can be generated using only the Clifford group which en¬ 
code the creation/annihilation operators of the molecular Hamiltonian and arbitrary single¬ 
qubit rotations which encode the electron energy (hjk) dependent time-evolution. For an 
explicit example of such a decomposition, see the Appendix. The coefficients in the molec¬ 
ular Hamiltonian can vary greatly from molecule to molecule and, at first, it seems that 
arbitrary rotations are required. We can remove this requirement by noting that arbitrary 
rotations can be generated from the repeated application of only two rotations. Almost 
any pair of rotations will work with the exception of two rotations that are symmetry el¬ 
ements of the same point group. The canonical choice in quantum information is H and 
T = exp(f7r/8)f?z(7r/4) also referred to as the n/S gate. The T-gate, which does not pre¬ 
serve the chiral octahedron, plus the Clifford gates is a common universal gate set. Note 
that T 4 = S 2 = Z. It has long been known that any rotation can be efficiently simulated 
using these gates m but, as we discuss in this review, recent work has drastically reduced 
the cost. 

Imagine a situation where qubit 1 can only have Clifford gates act on it (including 
the CNOT) and qubit 2 can have Clifford gates and rotations around the z-axis by an 
arbitrary angle R z {9). We can effectively perform R z {9) on the first qubit by using the 
important quantum primitive of gate teleportation pZ5J. The procedure, depicted in figure 
i>, works as follows: qubit 1 is in an arbitrary state. Qubit 2 is prepared in the state 
|+) = l/\/2(|0) + |1)) and the rotation R z {9) is applied. We will refer to this state as 
|A(6 1 )). Next CNOT(l,2) is applied and then qubit 2 is measured in the computational 
basis. The net effect is as follows: if qubit two is measured in state jO), then R z ( 9 ) was 
applied to qubit 1. If qubit two is measured in state 11), then R z (—0) was applied to qubit 
2. In this case, we can gamble and attempt to apply R z (29) and continue this procedure 
until success [26] . 

If 9 = 7t/4, we do not have to gamble. If we measure qubit 2 to be in the state 1, we 
have performed the undesired rotation R z (— vr/4). This can be corrected deterministically 
by applying the Clifford group operator R z (ir/2) = e~ m ^ 4 S to the first qubit (see figure 
r2p). We can then use our ability to perform R z {jt/A) to create a deterministic circuit for 
implementing R z (tt/8). This process can continue and it is often useful to consider gates 
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Figure 2: In this figure, we use the language of quantum circuits [18j to show the imple¬ 
mentation of general rotations, Rz{0 ), and the deterministic implementation of the rotation 
Rz( 7t/ 4) by the use of gate teleportation. Single black horizontal lines correspond to qubits. 
Horizontal double black lines to classical bits. Vertical lines represent control with a black 
(white) dot signifying control on the qubit or bit being 0 (1). Time flows from left to right, 
(a) Circuit for the implementation of the gate Rz(@)- Teleportation is attempted by the 
application of the rotation Rz(0) to the ancillary state |+), a CNOT to induce entanglement 
between the input state \ip) and the ancillary state |+), and the measurement of the ancillary 
state in the computational basis. If the measurement outcome is 0, then the desired gate 
was applied and the circuit terminates as indicated by the perforations of the quantum wire 
containing the input state. If the outcome of the measurement is instead 1, then the rotation 
was applied in the undesired direction about the axis and a correction rotation, Rz(26) , may 
be attempted in a similar manner. We may continue to apply these corrections until we get 
the desired measurement result, (b) Deterministic teleportation of the rotation, Rz( 7r/4), to 
the input register | ip) from the utilization of magic state, |H) = |H(7 t/ 4)). The first circuit 
consists of the unknown data state \i(j), the magic state |H) = |H(7 t/ 4)), a CNOT, a mea¬ 
surement in the z basis, and a measurement dependent Clifford operation, Rz( 7r/2). In the 
second circuit we write |H) in terms of the |+) and the gate Rz( 7r/4). In the third circuit, 
we have pushed the single-qubit gate through the CNOT, transforming it into a two-qubit 
unitary Rzz{ tt/ 4). In the fourth circuit, the CNOT can be removed since the identity and 
the NOT operator, A", have the same action on |+). In the fifth circuit, we note that the 
Z measurement commutes with the operator Rzz and transforms Rzz into a measurement 
controlled Rz- Finally, we observe that the measurement value does not change the overall 
operator applied to the first qubit. 
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that are of the form Rz( 7r/2 fc ) as well. 

For a physical qubit, for example a trapped ion or the nuclear spin of a molecule, the 
difference between rotations of different angles corresponds to simply a change in the pulse 
area. There is no greater difficulty in applying Rz( 7t/ 32) or Z. Many quantum systems have 
shown single-qubit gates with error probabilities on the order of ICC 4 or less per gate [2Thf33] . 
If these error rates can also be achieved with two-qubit gates, then a quantum algorithm with 
> 10 4 gates will not complete successfully, limiting the application of quantum computers to 
quantum chemistry. The solution to this problem is fault-tolerant quantum error correction. 


Quantum Error Correction and Preferred Universal Gate Sets 

For classical computers, error correction is possible due to the discrete nature of the states 
and the ability to copy information. At first glance, quantum computers and error correction 
seem mismatched since quantum states form a continuous space and the copy of quantum 
systems is forbidden by the no-cloning theorem [34]. The insight that allows for quantum 
error correction is that, instead of having discrete states, one can have discrete errors m- 
Each error shifts the data from one subspace of the system to another. A set of measurements 
can be applied that reveals which subspace is occupied without disturbing the data. The 
collapse of the state onto specific subspaces after measurement transforms the possible range 
of continuous errors into a set of discrete errors. 

For adequate partitioning into subspaces that can distinguish errors, a single logical qubit 
is encoded into multiple physical qubits. If the probability of an error on each physical gate is 
below a certain threshold, then arbitrarily long computations can be performed if sufficiently 
powerful error correction is used pTb! iHT| . The threshold theorem is a key result of quantum 
information and the value of the threshold depends strongly on the codes and error models 
but, for realistic noises and qubit geometries, it has been calculated to be as high as ~ ICC 2 
[381 [39] . Currently two physical systems, trapped-ion qubits and superconductor qubits 
have demonstrated a set of universal operations and measurements with error rates below 
this threshold [3IJ [40] . 

When considering a code, both the resources required for memory and for computation are 
pertinent. For every code, there is a small set of gates that can be implemented transversally. 
Transversal operations act on each physical qubit separately in the logical qubit. This implies 
that the resource overhead of transversal operations is equivalent to the cost of the encoded 
memory in terms of the numbers of physical qubits. A common example is where the logical 
operation is the same physical operation on all qubits. Some examples of quantum error- 
correcting codes and their transversal gate sets can be found in 0B. For many codes these 
correspond to the Clifford operations, for some codes only the Pauli operations, and for 
others only rotations of certain angles about the Z axis. For single-qubits, ignoring the 
global phase, these operations correspond respectively to the point groups O , D 2 , and D 2 k. 

For an important family of codes the Clifford gates are transversal and the costly gate is 
the T gate m This gate can be implemented by teleportation if the magic state |A(7t/4)) = 
T |+) can be prepared following figure^. At first glance this seems like an avoidance of the 
problem. Clearly to generate |A(7t/4)), we need to apply T. However, we can use a different 
code family to build a distillation circuit that allows us to efficiently generate high-quality 
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|A(7r/4)) states from multiple noisy copies. 

T Gate Distillation 

Bravyi and Kitaev developed the first method for T gate distillation based on the [15,1, 3] 
quantum Reed-Muller code to distill the magic state |A(7 t/ 4)) [19]. This code takes 15 noisy 
ancillary magic states, with error rate p, and encodes them into 1 qubit of information with 
a failure in the distillation circuit when 3 or more errors occur during the encoding. By 
performing the error check measurements and corrections in the quantum error correcting 
routine, one may iteratively project the noisy set of magic states onto one higher fidelity 
magic state, with an output error rate O (p 3 ), resulting in a noisy to distilled state ratio of 

i5-to-i m- 

The distillation circuit for this protocol is depicted in figure [3] Notice the circuit employs 
only Clifford operations except for the input magic states used for deterministic T gate 
teleportation as in Fig. Further rounds of distillation can be performed on distilled 
states to suppress error rates up to an arbitrary precision and, thus, allows for fault-tolerant 
implementation of the non-Clifford T gate. 

Meier, Eastin, and Knill (MEK) utilized a [10, 2, 2] error detection routine within the 
distillation circuit [43j. This code has a higher ratio of noisy to distilled states of 10-to-2 at 
the cost of a reduction in error suppression of the distilled states of 0(p 2 ). Dependent on 
the target error rate for the distilled state, the 10-to-2 technique can more efficiently produce 
distilled states than the 15-to-l protocol. There is a regime of low input error rates for which 
the quadratic suppression of errors from 10-to-2 technique is sufficient and extra resource 
cost for cubic error suppression via the Bravyi-Kitaev distillation is deemed unnecessary. For 
high-initial errors and low target errors, initial implementation of the 15-to-l protocol for 
maximum error suppression followed by the 10-to-2 protocol to minimize ancilla overhead 
can result in a favorable resource scaling relative to using a single protocol [43]. 

Bravyi and Haah introduced a method with a ratio of noisy to distilled states of 3k + 8 
to k for |t4(7t/ 4)) with output error rates comparable to that of Meier-Eastin-Knill, O (p 2 ) 
[44]. Furthermore, they introduced a systematic method of generating distillation protocols 
by realizing an equivalence between a set of quantum error correcting codes that admit a 
transversal T gate and a family of matrices known as triorthogonal matrices. The nonzero 
elements of these matrices in a row denote the set of simultaneous checks performed on 
the physical qubits per step of the error correction routine. Therefore, generating new 
distillation routines translates into generating new triorthogonal matrices. As with the 10- 
to-2 distillation, the utility of this technique is substantial at lower error rates and, therefore, 
admits serial use of distillation protocols due to constraints on the input error regime in which 
the resource gains of this technique outweighs the gain in output state precision [44] . 


Gate Decomposition Techniques and Minimizing Non-Clifford Rotations 

We will now present three general methods for generating fault-tolerant arbitrary unitary 
transformations. The first method uses a minimal, universal, elementary gate set (typically 
Clifford operations and T gates) and approximates arbitrary rotations through successive 
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Figure 3: The Bravyi-Kitaev 15-to-l magic state distillation protocol. The input to the 
circuit is 16 easy to prepare input states, |0) and |+), and 15 low-quality magic states with 
error p, and the output is 1 magic state with an output error rate of O ( p 3 ). The circuit is 
composed of Clifford operations and Pauli measurements. It implements 15 deterministic T 
gate teleportations (figure^). Subsequent rounds of distillation can be performed on the 
output states to suppress the errors of the magic state to arbitrary precision, thus allowing 
fault-tolerant implementation of the non-Clifford T gate. 






















































































application of gates from this set. These methods minimize the number of T gates to reduce 
the overhead due to distillation. The second method, known as complex distillation, incor¬ 
porates an overcomplcte gate set which is composed of Clifford operations and incremental 
rotations about the x, y , or z axis on the Bloch sphere to construct arbitrary rotations. The 
third method is sequential probabilistic application of gates through teleportation. Table [T| 
illustrates the decrease is resource costs from the progress in these techniques. 

Optimization with a Minimal Gate Set 

Gate decomposition techniques critically rely on a fundamental result in quantum computing 
known as the Solovay-Kitaev theorem. Given a set of single-qubit gates (e.g., if, T, and 
T t) that generate a “dense” set of unitary gates, it is possible to approximate any unitary 
gate to an arbitrary precision (5), by an efficient number of operations, oc log c (l/<5) where 
c is a constant m- The challenge is how to compile these elementary gates into accurate 
approximations of arbitrary unitary operations. 

We will now outline the progress in methods of approximating arbitrary unitary rotations 
from a set of Clifford and T gates. The Solovay-Kitaev algorithm presents a step-wise 
procedure for generating arbitrary gates B5j. First, a library is generated of all rotations up 
to a sequence length of elementary gates l. The target rotation is compared to the closest 
rotation in the library. The difference between these rotations is a small rotation that can be 
generated from the library taking advantage of the non-commutivity of rotations. This can 
be done in a recursive manner to generate arbitrarily accurate gates without ancillary qubits. 
Further optimization methods have been implemented to the Solovay-Kitaev algorithm with 
improvements that reduce the level of recursion required for the calculation of the sequences 
|46j. The final elementary gate count also depends critically on the size of the library as 
demonstrated by Bocharov and Svore m- 

The best gate counts will occur for infinitely large libraries. However, this approach be¬ 
comes impractical for high-precision approximations. For moderate precision, Fowler found 
optimal sequences by a “brute force” space search that utilizes information about previously 
applied gate sequences to selectively screen for the next gate sequence [46j @8]. The precision 
at which this method can approximate gates is limited as the space search is exponential 
in gate sequence length, but this result was nevertheless incredibly insightful as it illumi¬ 
nated the significant gap in gate sequence length between optimal decompositions and the 
Solovay-Kitaev result. 

An alternate approach for developing an optimal decomposition required the use of an¬ 
cillary states that could offer small corrections to the approximated state through a method 
known as phase kickback [36] • From these methods, a representation developed where the 
problem is mapped to a unit field and optimal gate sequences that can be obtained without 
the use of ancillary states [39lf52| . Non-deterministic methods of approximating such rota¬ 
tions have also been developed which make use of a minimal set of gates and ancillary states 
which are measured to iteratively project the input state to the required state within a given 
precision [211 loIH 166] . 

Significant progress has recently been sparked by work performed by Kluichnikov, Maslov, 
and Mosca (KMM) [56]. This work presented a rigorous proof that any single-qubit unitary 
operator can be exactly decomposed into a set of single-qubit Clifford and T-gates, with 
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one |0) ancilla state, if and only if the matrix elements of the unitary operator belong to a 
special algebraic group and fulfill constraints on the assignment of the matrix elements to 
guarantee unitarity. Furthermore, it was shown that this decomposition was efficient. Giles 
and Selinger extended this proof to include n-qubit unitary transformations and eliminated 
the need for ancillary qubits in the exact decompositions 153 This result is powerful as 
it added intuition into the search space of potential circuits resulting in a deviation from 
random circuit search methods. 

The question remained: given a desired rotation within a given precision, how does one 
find a rotation with an exact decomposition that lies within that precision? Furthermore, 
can one find an exact unitary decomposition, within that precision, which is minimal in 
the use of the costly T gate? Various algorithms for searching the space of unitaries with 
exact decompositions have been implemented for the generation of these circuits [^9H52] . 
These methods differ from the search methods implemented by Fowler in that, instead of 
searching over gate sequences, they search for the entries in the single-qubit unitary that 
simultaneously satisfy the norm equation to guarantee unitarity and are members of an 
algebraic ring, the cyclotomic integers, to ensure the existence of an exact decomposition. 
Initial attempts at these optimizations resulted in decompositions that did [59] and did not 
[50] require the use of a few ancillary |0) states in the circuit. Additional algorithms resulted 
in an optimal solution in the number of T gates with ancilla [51], and Ross and Selinger 
recently showed optimal solutions in gate decompositions that are exact and require no 
ancillary qubits m Ironically, the efficiency of the decomposition depends on the ability 
to factor numbers, a task at which quantum computers excel [58]. Fortunately, probabilistic 
versions of this search algorithm can be implemented in polynomial classical runtime and 
can achieve the third-to-optimal decomposition sequence, with T gate counts comparable to 
the optimal solution [52] - 

The Utility of an Expanded Gate Set from Complex State Distillation 

Complex distillation protocols [59ff6T] utilize an overcomplete elementary gate set consisting 
of Clifford operations and a set of rotational states about a given axis. These rotations 
about an axis are produced by a distillation protocol that takes magic states, such as those 
constructed from the routines above, as an input and generates a “harder” rotation. While 
employing subsequent distillation procedures on distilled states appears to lower efficiency, 
the rationale behind these techniques is founded in the fact that the minimization of the 
number of ancilla required per magic state does not necessarily imply more cost-efficient 
computation as the above protocols are ignorant to the compiling cost of the computation; 
an issue addressed by these complex distillation techniques. Indeed, there is gate number 
threshold above which generating a “difficult” rotation directly may benefit over successive 
applications of a set of cheaper rotations when approximating an arbitrary rotation. As an 
example of how this may reduce a gate sequence length for an arbitrary desired rotation, if 
one has access to the rotations Rz( 7r/2 fc ), one can approximate any arbitrary Rz(0) by digital 
compilation. This can be especially useful for algorithms that require on-the-fly compilation 
of unitary gates based on measurement of qubits [61, E2J. 

These protocols differ in the expanded gate set used and the distillation used to gener¬ 
ate the expanded magic states. What they do hold in common is that the expanded set 


10 


of rotations about an axis are generated recursively, ^ R{9 \) —> R{6 2 ) —> ... —> R z (0 k )j, 
with finer rotations requiring more distillation steps for all techniques and these rotations 
are applied to data via gate teleportation. 

Duclos-Cianci and Svore [59] use an elementary gate set that includes Clifford gates and 
Ry rotations. The resource states used in the teleportation were constructed from distillation 
of | Hi) magic states, | Hi) = Ry(2Qj) |0) where 0* is defined by cot0j = cot l+ 1 7 r/ 8 . Note that 
| Hq) = |+). An additional routine utilizing a set of additional rotations {|-0°), \i'pl) , |-0?)} 
generated Clifford circuits with |iJj)’s as an input creates a larger, denser set of rotations 
partitioning the xz plane. The work presented a systematic method of generating dense 
states about a plane on the Bloch sphere. The rotations generated by the magic states are 
not uniformly partitioned and although they can be combined to generate any angle, there 
is no natural decomposition. 

Landahl and Cesare modified the previous approach to complex distillation techniques in 
two ways: first they addressed the non-uniform partitioning of the yz plane by performing 
Rz ( 7 r/ 2 fc ) rotations and also implemented a distillation routine with a “top down” recursion 
[60]. In the “top down” recursion method, additional R z ('rr/2 k ~ 1 ) input states are needed 
if the R z (vr/ 2 fc ) teleportation is faulty and even less R z (vr/ 2 fc ~ 2 ) states are required as 
they are applied upon consecutive failures of the previous two teleportation attempts which 
probabilistically reduces the number of traditionally distilled states required for gate tele¬ 
portation because the teleportation must fail k — 3 times before a T |+) magic state must 
be prepared ED]- The technique used shortened Reed-Mullcr codes that require input 
magic states of increasing quality to achieve the smallest rotations. 

Duclos-Cianci and Poulin overcame the limits of using Reed-Muller codes by the im¬ 
plementation of a modified 10-to-2 distillation circuit [?3] to generate magic states 11*,) = 
ZSH |A(27r/2 fc )} for Ry in/2 k \ rotations. Conditional on the measurement, the distillatoin 
procedure transforms two noisy \Y k ) states into two quadratically improved \Y k ) states. The 
method requires Clifford operations and the application of a parity measurement, which 
consumes 16 magic states of the type I 3 ) and 1 magic state of the type \Y k _\) j 6 Tj . The 
similarity of the distillation for all k allows the procedure to work with magic states of a fixed 
initial accuracy regardless of k. For an input accuracy of 1%, 10), a Clifford state, serves 
as a sufficiently accurate input state. The authors find a slightly lower resource overhead 
in terms of non-Clifford input states when substituting |0) for input \Y k ) when k > 8 than 
for k > 3. This method has comparable resource costs relative to T-gate only compilation 
methods for arbitrary rotations, but substantial savings for rotations by angles 2n/2 k . 

Non-Deterministic Application of Rotations 

The third set of methods can be thought as random-walk with correction teleportation tech¬ 
niques. Teleportation is attempted with ancillary states that are measured to detect the 
success of the application of the gate, similar to the teleportation scheme shown in figure 
[2p. If the gate fails, recursive corrections, such as those we have seen from complex dis¬ 
tillation techniques, can be applied until the desired gate is applied. An early example of 
this technique was shown by Jones et al. which implemented programmable ancilla rota¬ 
tions (PARS) [2l|. This method required non-Clifford ancillary states that were originally 
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expensive to generate but distillation techniques have made the production of Fourier states 
efficient [63]. Recent methods have eliminated the need for these complex ancillary states to 
non-deterministically achieve arbitrary rotations [S3] EHJ • The most recent non-deterministic 
technique by Bacharov et al. utilizes the knowledge of previous rotation attempts in a man¬ 
ner that probabilistically guarantees the termination of the circuit after a finite set of gates 

H- 


Method 

T Gate Count for Precision 5 

Non-Clifford States 

Solovay-Kitaev a [J5] 

s, (log’-'-" (1/<S)) 

86 

8.0 x 10 5 

Fowler Search a [12 08] 

2.95 log 2 (l/<5) + 3.75 

1.3 x 10 2 

3.4 x 10 2t 

PARs 5 [2Tj 

s 2 (log (l/d)) 

6.0 x 10 2 

1.4 x 10 3tt 

KMM“ [ST] 

3.0671og 10 (1/(5) — 4.322 

88 

8.2 x 10 2 +tt 

Floating-Point 5 [53] 

8 log 2 (1/(5) + 1.14 log 2 (10 7 ) 

> 2.9 x 10 2 

> 2.9 x 10 3 

Ladder States c [59] 

- 

1.1 x 10 2 

4.3 x 10 3 

CISC C K3) 

- 

7 

34 

RS a [52] 

3 log 2 (1/5) + log 2 (log 10 (1/(5)) 

1.1 x 10 2 

1.1 x 10 3 

rus 5 [2SIEU 

1.15 log 2 (1/(5) 

42 

4.2 x 10 2 

dcp c pm 

- 

7 

34 

PQF b PH] 

log 2 (1/5) + log 10 (log 10 (1/5)) 

38 

3.7 x 10 2 


Table 1: Resource states required for non-Clifford gates for the various methods discussed 
in the review. Techniques are presented approximately chronologically. Two assessments 
of the number of non-Clifford states are presented: a target gate precision (<5) of ICC 2 in 
the left column and 5 = ICC 20 on the right except for the Fowler (f), PARs (ff), and 
KMM (fft) which are evaluated at a precision of 10 -6 , 10 -5 , and ICR 15 , respectively. The 
lower bounds on the floating-point method are given as the T gate assessment target is 
angle dependent {6 = a x 10 -7 ). The superscripts denote the category of the technique as 
presented in the text: (a) Clifford and T decompositions, ( b ) non-determinisitc methods, 
(c) complex distillation. The estimates above assume the availability of sufficiently precise 
magic states. This assumption increases the advantage of complex distillation techniques. 
The Ladder State value relative to CISC and DCP is inflated as the number includes the 
error-free distillation cost of the higher-order states. The advantage of DCP over CISC is in 
the distillation procedures for noisy magic states (see text). The scaling factors Si and s 2 
depend on the specific algorithm but are > 1. 


FUTURE DIRECTIONS 

Without functioning devices, adequate estimation of resource costs of quantum simulation 
remains difficult. When one considers the wide use of heuristics and approximation in 
quantum chemistry by classical computation, it seems likely that a similar explosion of 
methods will occur when a real device is present. 

The advantage of a quantum chemical calculation on a quantum computer is the low cost 
of expanding configuration space. The disadvantage is that, for the canonical algorithm, 
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the energies are calculated by solving the dynamics problem. This suggests that quantum 
computers may have a larger advantage for calculating dynamic properties [64j and also that 
new algorithmic ideas, such as variational solvers [63], may lead to better ways to retrieve 
energies. 

We have outlined resource efficient methods of approximating arbitrary qubit rotations 
that will be required in computing the elements of the Trotter decomposition of the time 
evolution of a molecular Hamiltonian. An open question is how to decompose, within a given 
precision, the continuous time evolution of the Hamiltonian into a set of Trotter steps in a 
manner which is efficient (logarithmic or polylogaritlnnic) in resources and time during the 
simulation. Algorithms have been developed which have displayed such decompositions and 
exploit trade-offs between the number of qubits and the number of parallel computational 
steps [66Tf69j . Connecting the Trotter compilation to rotation composition and distillation 
could yield further reductions in the total number of gates, a notion partially investigated by 
Hastings et ah im Furthermore, an understanding of the required precision of the rotations 
implemented in the quantum chemistry algorithms to achieve meaningful results is critical 
for implementation into physical devices. 
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APPENDIX 

In this section we will show how to map the dynamics of the molecular Hamiltonian (Equa¬ 
tion [Tj) , more specifically an element of the Trotter decomposition, onto a set of Clifford 

operations and arbitrary qubit rotations. Given a Trotter step exp we can uti¬ 

lize the following identity: Let g be the Clifford operator such that gZjg^P, where P is a 
Hermitian Pauli operator (ZiZjXkXi for example), then the unitary Rp(9 ) = exp(— i6P/2) = 
g exp(—idZi/2)g^ = gR z (0)g\ where R z (0) is the rotation of qubit j about the z-axis by 
an angle 9. The key result is that the parameters of the Hamiltonian are encoded into the 
simulation by single-qubit rotations. 
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We no show how to find the Clifford operator g that transforms a rotation into a piece of 
the molecular Hamiltonian. Consider the term hj k {&^c k -\-c k Cj) in the molecular Hamiltonian 

c]c k + c f k £j 


where hj k is real. We wish to implement the unitary operator exp y—i(hj k t/K) 
for some time t. Utilizing the Jordan-Wigner transformation, we have 


N 


st _ 


= (Xj - iYj) n X,- 


1 = 3 +1 


The long chain of Z operators guarantees the correct anticommutation behavior. Without 
loss of generality, assume k > j then 


c)c k = (X 3 - iYj) ( J] zA (X k + iY k ) 
\i=j+1 / 


and we see that 


k -i 


fc-i 


+ eldj = Xj I H z,\ x t + yA n Z‘) U. 


\l=j+1 


\l=j +1 


is the sum of two Pauli operators that commute. 

The next step is to build the Clifford operators to transform a single-qubit Z operator into 
these (k — j + l)-body spin operators. We use three identities: CNOT(/, m)Z m CNOT(/, m) = 
ZiZ m allows us to grow a chain of Z operators, HZH = X transforms Z’s to X’s, and 
SHZHS t = Y transforms Z’s to Y' s. 

We apply these identities to the problem of finding a Clifford operator gxx that satisfies 
the equation gxxZ k (/ X x = Xj ^n?=/+i X k . The solution is to grow a chain of Z’s from 

k to j, (nt; CNOT(Z, l + 1) j Z k (nf - 1 CNOT(/, l + 1) j = nf=y an d then apply two 

H' s; gxx = HjH k CNOT(/,Z + l)j . A similar analysis for the other Pauli operator 

reveals g Y y = SjS k g X x- 

Combining all of the pieces we have 


exp 


-i(h jk t/h) 


st 


t; 


dc k + c' k Cj 


gxxRz k {2hj k t/h)gl cx g YY Rz k ( 2hj k t / h)g\r Y 


which is a unitary operator made from Clifford operations and a rotation about the z axis 
set by the molecular Hamiltonian and the time of simulation. Although the details of the 
transformation differ, this gate set is also all that is required for the Bravyi-Kitaev mapping 

HSj. 
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